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Abstract 

We study Brownian motion driven with both conservative and nonconservative external forces. 
By using the thermodynamic approach of the theory of Brownian motion we obtain the Fokker- 
Planck equation and derive expressions for the Fluctuation Theorem in local equilibrium and in 
quasi-equilibrium. In local equilibrium the expressions we obtain coincide with previous results. 
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I. INTRODUCTION 

Since the first formulation of a Fluctuation Theorem (FT) llfl there has been an increasing 
interest in this subject. Further contributions [2] have broadened the scope of the physical 
situations where a FT arises. An extensive account of the literature can be found in Refs. 
[3] and 4]. The search for certain symmetries or rules in the different systems or situations 
we can find far from equilibrium is the major element of interest in the different expressions 
of the FT. It seems that the FT play a role similar to the partition function in equilibrium 
statistical mechanics. 

These theorems are concerned with the relation between nonequilibrium measurements 
and certain thermodynamic magnitudes such as entropy, heat and free energies, 4j. They 
involve pairs of nonequilibrium trajectories in the phase space, one and its reverse, and their 
corresponding probabilities. 

Our contention here is to clarify the thermodynamical roots of the FT, determining their 
scope for systems with Fokker-Planck dynamics. We will obtain the expression of the FT in 
local equilibrium and in states of quasi-equilibrium. 

The paper is organized as follows. In section 2, we perform the thermodynamic analysis 
of an ensemble of non-interacting Brownian particles under the action of both conservative 
and nonconservative forces, deriving the Fokker-Planck equation. In section 3, we derive the 
configurational relaxation equations when the system is in local equilibrium and in quasi- 
equilibrium and obtain the parameters which characterize these states. Section 4 is devoted 
to the derivation of the FTs. Finally in section 5 we formulate our main conclusions. 



II. THERMODYNAMIC ANALYSIS 



Consider an ensemble of non-interacting Brownian particles subjected to a potential field 
V(x) which is initially in equilibrium with a heat bath at temperature T Q . In the canonical 
ensemble, the system is distributed according to 

(1) 

where H {V) = (l/2)mu 2 + mV(x) is the Hamiltonian and T = (x,u) represents a point of 
the one-particle phase space spanned by the position x and the velocity u of the particle. 



p eq . ~ exp 



g(T) 
k B T a 
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Let us assume that at certain moment (t = 0), a nonconservative force f(t) is applied on 
the system. Then, the system evolves in time and can be described by using the nonequi- 
librium probability density p(T, t) and the nonequilibrium entropy of the system which is 
given through the Gibbs entropy postulate {^{(j] 

S(t) = -k B [ p(T, t) In ?^-dr + S eq .. (2) 

J Peq. 

where S eg . is the equilibrium entropy of the Brownian particles plus the bath. The variations 
in the probability density p(T, t), imply changes in the nonequilibrium entropy which can 
be obtained from Eq. (J2J) 

5S = -T I Kr,t)6p(r,t)dr. (3) 

The thermodynamically conjugated nonequilibrium chemical potential p(T, t) of the density 
p(T,t) has been defined by 

p{T,t) = k B T \n^^ + p eq ., (4) 

Peq. 

where p eq is the equilibrium chemical potential. Equation ([3]) is similar to the Gibbs' equa- 
tion of thermodynamics 

After the force f(t) has been applied, the variation in time of p(T, t) is given by the 
generalized Liouville equation 

J^(r,t) + v r (r,i)-Vr/>(r,t) 

where Vr(T,t) = (u, — VV r (x)) is the phase space velocity corresponding to the Hamiltonian 
flow, Vr = (V, d/du) with V the spatial derivative and J(T, t) constitutes a diffusion current 
in phase space. By using Eqs. ([3]) and (jSJ) and performing partial integrations assuming 
that the currents vanish at boundaries, one obtains that the total rate of change of the 
nonequilibrium entropy (j2J) is 



f = -£w>«>-£/jCr,«>!rtr,«)«ir, (6) 



where the quantity (f(t)u) = J f(t)up(T,t)dT constitutes the power, dw/dt, supplied by the 
external force /(£) which is dissipated as heat into the bath. Thus, assuming that [3] 

p(T,t) = 8(x(t)-x)8(x(t)-u), (7) 



we can write 

^ = /(*)*(*), (8) 

where dw is the amount of work done on the system in a time dt by the external force /. 

The second term on the right-hand side of Eq. ([6]), represents the change of the entropy 
due to a gradient of the chemical potential and therefore constitutes the entropy production 
a due to diffusion in the phase space 



a 



1 f d 

-J J(T,t)—fi(r,t)dT. (9) 



This quantity accounts for the internal dissipative processes and it is assumed, according to 
the second law of thermodynamics, that it does not take negative values: a > 0, [6j. 

Hence, by introducing the exchange of entropy due to the interactions with the environ- 
ment as 

d ext S 1 dw 

~JT = T ~dl 1 ( ' 

which has an undefined sign, we finally obtain 



dS d ex iS 
dt dt 



+ <j (11) 



which expresses the entropy balance between the exchange of entropy with the bath and 
the entropy generated in the irreversible processes established in the system when removed 
from its equilibrium state. 

To completely characterize the stochastic dynamics of the system, we must find the 
expression of the current J(T,t). According to Eq. ([9]) and following the rules of nonequi- 
librium thermodynamics {(|, this can be achieved by establishing a linear relation between 
the current J(T,t) and their conjugated thermodynamic force dfi/du. This relation can be 
expressed by 

J ™ = 4 1^> < 12 > 

where L is a phenomenological coefficient. By using the expression of the chemical potential 
given through Eq. (jlj) in (fl2l) . one obtains 

J ( r ' t » = -' '(^ + u )'<™ (13) 
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where we have identified mL/ pT a = 7, with 7 being the friction coefficient per unit mass of 
the Brownian particles. By substituting now Eq. ffT3j) into Eq. ([5]) we obtain 



_ p = _ Vup + _ [W (x) - /(f)] p 

d fk B T d \ 
+ — -7T + U )p (14) 



du \ m du 

which is the Fokker-Planck (or Klein-Kramers) equation governing the time evolution of the 
nonequilibrium probability density p(T, t) in the presence of a conservative potential V(x) 
and a non-conservative force f(t). 

III. LOCAL EQUILIBRIUM AND QUASIEQUILIBRIUM 

When the dynamics of the system can be characterized by two or more time scales, its 
long-time behavior can manifest different dynamical regimes depending on the existence (or 
not) of a local equilibrium state. Here, we will show the conditions that must be satisfied by 
the system in order to reach a local equilibrium state and how, when these conditions are 
not satisfied, the system enters in a quasiequilibrium state leading to a different long-time 
dynamical behavior. 

Assuming the presence of two time scales in which u is the fast variable, the dynamical 
processes associated to configurational changes in the system are related to the slow variable 
x. Thus, the long-time behavior of the system can be more conveniently described by writing 



p(T,t) = (p x (u,t)n(x,t), (15) 

where (p x (u,t), is a conditional probability density and n(x,t) = f p(T,t)du is a reduced 
probability density in x-space which evolves according to 

d f 

m—n(x,t) = —'V / mup(T,t)du. (16) 

This equation can be obtained by integration over u of Eq. f[T4"]) and thus implicitly defines 
the current J(x,t) = f mup(Y ,i)du. 

i) Local equilibrium - In order to achieve a local equilibrium state, the system must satisfy 
that <p x (u,t) ~ exp (— mM 2 /2/c j eT ), i.e. that the distribution of velocities is given by the 
equilibrium Maxwellian. In this case, by multiplying Eq. f|T4l) by mu and performing a 
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partial integration over w-space, we obtain the evolution equation for J(x, t) 

d 

— J(x, t) + 7 J(x, t) = -mn(x, t) [VV(x) - f(t)] 

-k B T Vn{x,t). (17) 

For times t ^> 7 _1 Eq. (JTTJ) gives 

J(x, t) ~ -7" 1 {mn(x, t) [W(z) - /(£)] + k B T Vn(x, £)} . (18) 

In order that our results be valid, here and henceforth we will assume that the time scale 
over which f(t) varies should be larger than 7 -1 . Otherwise, the time derivative of Eq. 
(ITTjl must be taken into account. Once we substitute Eq. (fl8"|) into Eq. (|T6|) we obtain the 
Smoluchowski equation in the presence of external forces 

-n(x,t) = DV l— n (x,t) [VV(x) - f(t)] + Vn(x,t)\ (19) 

with D = ksTo/m^ being the diffusion coefficient. The Smoluchowski equation f|T9|) admits 
a local equilibrium solution given by 



ni.eq.ip; f) ~ exp 



m 



[VV{x') - f(t)} dx' 



k B T 

Hence, the nonequilibrium entropy given through Eq. ([2]) reduces to 



(20) 



/Til X t) 
n(x,t)ln^^ L dx + S Leq Xt), (21) 
ni.eq. 

where Si. eq .(t) = —k B J n(x, t) ln[nz. eg .(x; f)/n eq ]dx + S eq . is the local equilibrium entropy. It 
is important to emphasize that through Eq. ( |20l) . the local equilibrium state is characterized 
by a probability density containing the temperature of the bath T . 

ii) Quasi- equilibrium - Released from the restrictive condition of local equilibrium, instead 
of obtaining the expression (fTBl) for the diffusion current J(x,t), at long time we find 

J(x, t) = -7- 1 {mn(x, t) [VV(x) - f(t)} + k B Vn(x, t)T(x, £)} , (22) 

where the local temperature of the system T(x, t) has been defined through a generalization 

nn 

of the equipartition theorem as the second moment of 4> x (u,t) [10], [11] 

k B T(x,t)=m J u 2 (fi x (u,t)du. (23) 
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The current given in Eq. (1221) can be rewritten as 
J(x,t) = -D(x,t) 



V/;(,-./) + 7 m ^ X J^ [V$(x,t)-/(t)] 



(24) 



k B T(x,t) 

where $(x, t) = V(x) + ksT(x,t) is an effective potential and D(x,t) = fcsT(x, t)m _1 7~ 1 
is the bare effective diffusion coefficient. Hence, by substituting Eq. fl24l) into Eq. (fT6l) we 
obtain 



5_ 

5f 



n(x,t) = V <^ D(x,t) 



The diffusion equation ( |25i) given above admits a quasi-equilibrium solution 



n qe ,(xj) ~ exp 



[V$(x',/)-/(t)]^ 



(25) 



(26) 



k B T(x'J) 

obtained when the probability current instantaneously vanishes. 

An estimation of the temperature T independent of the position can be obtained by first 
deriving the evolution equation of the temperature field T(x,t). Thus, by multiplying Eq. 
( TT4"j) by mu 2 and integrating in u one has 



1 d 

-/cs— n(x, t)T(x, t) = — Vn(i, t)h(x, t) — 
[W(z) - /(*)] J{x, t) - jk B n{x, t) [T(x, t) - T ] 



(27) 



where we have defined the heat flow as 



h(x,t) = -m I u 3 (f) x (u,t)du. 



(28) 



For times t ^> 7 1 , Eq. fl28|) reduces to 



ksn^x, t) [T(x, t) — TJ = —7 1 Vn(x, t)h(x, t)- 
7 -i [W(x)-/(t)] J(x,t). 



(29) 



In the particular case of local equilibrium J(x, t) = h(x, t) = which would lead to T(x, t) = 
T Q . If, on the other hand, we subsitute the expression of J(x, t) given through eq. ( 1221) into 
Eq. fl29l) and using Eq. fl29|) . up to order 7~ 2 we obtain 



(30) 



kBn(x, t) [T(x, t) — T Q ] = —7 1 Vn(i, t)h(x, t)+ 
7 " 2 {mn(x, t) [W(x) - f(t)f - k B T [W(x) - f{t)\ Vn(x, t) } 
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which after performing an average over x-space leads to 

k B T(t) ~ k B T + m 7 - 2 (W(x) - f(t)} 2 



(31) 



Here T(t) = (T(x,t)) and we have neglected a term k B T Q (V \yV(x) — fit)}). Relation (l3"Tj) 
means that for sufficiently large gradients and forces applied on the system, its temperature 
will in general differ from that of the heat bath. On the contrary, it must be emphasized 
that for small gradients and forces the quasi-equilibrium temperature T(t) given by Eq. 
( 13"Tj) reduces to the bath temperature T Q implying that the system reaches local equilibrium. 
This is precisely the hypothesis of small forces and gradients which is usually behind the 
local equilibrium hypothesis 12]. A correction to the temperature similar to ( 13TT) has been 
previously obtained in a different context in Ref. [lo| . 



IV. FLUCTUATION THEOREMS FOR LOCAL EQUILIBRIUM AND 
QUASIEQUILIBRIUM 

The transition from the intial state x Q to the final state Xt, can be modeled by a set 
of N coupled unimolecular chemical reactions having x as their reaction coordinate. This 
approach is based on the fact that the kinetic or rate equation corresponding to a uni- 
molecular chemical reaction represents a gain and loss process which can be interpreted 
probabilistically as the result of the balance of two opposite probability currents. 

To compute the ratio between the forward and reverse probabilities of a path between 

x and x t , we define a partition ti,t 2 , , tjv+i of the entire time interval [0, t), with ti — 

and ijv+i — t, which divides the process in N steps. The initial state for these reactions 
coincides with x Q and the final state with x t , while the intermediate states correspond with 
Xt 2 , Xt 3 , .... Thus, given the probability n, of being in the state Xt % at time ti, the elementary 
kinetics is given by the set of equations 

_l n . = J*-J < - 1 ; i = l, ,N + 1, (32) 

dr 

where 

J 4 = v* hi n i+1 - v[ ii+1 m; J° = J N+1 = 0. (33) 

Here, the symbols F and R stand for the forward and reverse transitions. On the other 
hand, it must be stressed that the values of the rate constants v^- r {t) depend on whether 
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the system is in a local equilibrium or a quasi-equilibrium state. Therefore, since the reaction 
constant j-jj characterizing the asymmetry of the reversible reaction is defined by the ratio 
between the rate constants, this will depend also on the state of the system. According to its 
definition, r^- gives us the ratio of the probabilities of the forward and backward reactions. 

In the quasi-stationary state J 1 = J 2 = = J N = and from Eqs. (1321) and (1331) we 

find 



lim 



71; 



nUi(f) 



J i,i+1 



(/) 



r 



i,i+l 



(f) 



(34) 



— n*(r) ~ n?(/) " t& li4 (/) 
where the upper index o stands for the long time value of n^, and r iti+ i(f) is the partial equi- 
librium constant corresponding to the i-th reaction (step), which depends on time through 
/. Thus, for the forward and reverse driven processes between t\ and ijv+i, we obtain the 
general relations 



ri, N +i(f) = ri i2 {f)r 2 ,3(f)...r N)N+1 (f) 



n° 2 (f)n° 3 (f) n° N+1 (f) 



TjN n i+l 
il i=l 



(/) 



(35) 



</)</) n%(f) i=i <(/) 
where now n° will have different dependencies on the temperature depending on whether 
the system is in a local equilibrium or a quasi-equilibrium state. 

i) Local equilibrium. By making use repeatedly of Eq. ( 1201) in Eq. (1351) one obtains 
that the ratio between the forward prob F (path) and reverse prob R (path) path probabilities 
is given by 

prob F (path) 



prob R (path) 



n,N+i(f) 



exp 



exp 

AF 

k B T Q 



m 

ksTo Jo 



-VV(y) + f(r)]ydr 



exp 



W 



1 



kEiTo 



(36) 



where Wf = m J Q * f(r)ydr is the work due to the nonconservative force in the forward 
path. Eq. (|36|) constitutes the ratio between the forward and backward path probabilities 



providing us the Fluctuation-Theorem at local equilibrium 
(IE]) and ( fTOl) . we can rewrite Eq. (|36l) in the form 

prob F (path) 



131 ] . Additionally, by using Eqs. 



e -AF/k B T 0e &/k B 



(37) 



prob R (path) 

where a = J Q T dc ^ s dt is the total entropy produced by the non-conservative force along the 



process. Eq. ( 1371) expresses that the ratio of forward and backward probabilities is related 
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to the entropy production [14| . Nonetheless, it should be pointed out that a is the entropy 
change of the system with its surroundings and therefore is of undefined sign. Another 
consequence of Eq. ( 1551) is the following: 



exp 



W 



k B T Q 



'y^ j prob F (path) exp 

path 

-AF 



W 



f 



k B T Q 



s ^j>rob R (path) exp 

path 

~—AF~ 
= exp — — 

k B T Q 



k B T a 



(38) 



which constitutes the so-called nonequilibrium work relation [15f. Despite of the common 

nn 

belief, Eq. ( 158"]) is only valid at local equilibrium as was stated previously [16|, [17] • 

It is convenient to emphasize that the Fluctuation Theorem adopts so simple and elegant 

expressions like Eqs. f[3"oT) and f|3T|) because we have derived them in the conditions of local 

equilibrium, i.e., when the fluctuation-dissipation regime is satisfied. Otherwise more fuzzy 

expression is obtained, as we will show next. 

ii) Quasi- equilibrium. Following a similar analysis as in the previous case, from the quasi- 

equilibrium solution given in Eq. (126]) we find 

prob F (path) 



prob R (path) 



r l,N+l(f) 



exp 



m 



[-VV{y) + f{r)]ydr 



(39) 



Uo k B T(r) 

giving us the ratio between the forward and backward paths probabilities at quasi- 
equilibrium. Eq. ( |39l) constitutes a generalization of the Fluctuation Theorem to the case in 
which the system not yet relaxed to local equilibrium, but it is still at a quasi-equilibrium. As 
we have mentioned previously, in far from equilibrium conditions, it is not generally correct 
to assume that temperature of the system is that of the heat bath, and therefore in far from 
equilibrium conditions Eq. (139]) must be used in place of (136]) . From Eq. (139]) it is not pos- 
sible to obtain a relation similar to (138]) because of the dependence of the quasi-equilibrium 
temperature in / and time along the path. 



V. CONCLUSIONS 

In this paper based on the thermodynamic theory of Brownian motion Q], we have 
derived the Fokker-Planck equation for a system of Brownian particles subjected to both 
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conservative and nonconservative forces. In this scenario we have obtained expressions for 
the FT in local equilibrium and in quasi-equilibrium. Local equilibrium coincides with the 
fluctuation-dissipation regime where the fluctuation-dissipation theorem is satisfied and the 
intensive parameters characterizing the state of the system coincide with those of the bath. 
In this case our expressions for the FT agree with previous results in the literature. 

Nonetheless, when arbitrary forces remove the system from equilibrium with the bath and 
keep it out of equilibrium modifying the intensive parameters characterizing the state of the 
system, the fluctuation-dissipation theorem is no longer valid. However, if these external 
forces vary quite slowly in time, one can define a state of quasi-equilibrium. Fluctuations 
around this quasi-equilibrium state are characterized by a generalization of the FT in which 
the existence of quasi-equilibrium state is taken into account. The theory we have presented 
here enables us to derive steady state fluctuation theorems for arbitrary forcing, something 
which we will do in a latter work. 
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